October 1, 2012 



NNLL Resummation for Jet Broadening 



Thomas Becher and Guido Bell 

Albert Einstein Center for Fundamental Physics 
Institut fiir Theoretische Physik, Universitat Bern 
Sidlerstrasse 5, CH-3012 Bern, Switzerland 



Abstract 

The resummation for the event-shape variable jet broadening is extended to next-to- 
next-to-leading logarithmic accuracy by computing the relevant jet and soft functions at 
one-loop order and the collinear anomaly to two-loop accuracy. The anomaly coefficient 
is extracted from the soft function and expressed in terms of polylogarithmic as well 
as elliptic functions. With our results, the uncertainty on jet-broadening distributions 



is reduced significantly, which should allow for a precise determination of the strong 



D , coupling constant from the existing experimental data and provide a consistency check 



on the extraction of a s from higher- log resummations of thrust. 
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1 Introduction 



Event-shape variables characterize geometrical properties of final states observed in high- 
energy collisions. They are the simplest class of observables which probe properties beyond 
the total cross section and have been measured with high accuracy at LEP and other e + e~ 
machines. Due to their inclusive nature, they have little sensitivity to hadronization effects 
and can be computed perturbatively. The canonical event shape is thrust T. To evaluate the 
thrust of an event, one first determines the thrust axis Ht, which is the direction of maximum 
momentum flow. The thrust is then defined as the fraction of particle momentum flowing 
along the thrust axis. The variable we consider here is the jet broadening, which measures the 
momentum flow orthogonal to the thrust axis. The total jet broadening is defined as [Tjr2] 



By using the thrust axis to split the event into left and right hemispheres, one can define left 
and right broadenings bi and 6r. The total broadening is the sum of these, bx = b^ + 6r, 
and the wide broadening is defined as the maximum of the two, by/ = max(6^, bn). One often 
normalizes the broadening to the total momentum flow which, for massless particles, is equal 
to the center-of-mass energy Q of the collision. We use capital letters to denote the normalized 
values, for example Bt = br/Q- 

For events consisting of two low-mass jets recoiling against each other, the thrust is near 
the maximum value T — 1, while the broadening is very small. In this region, higher-order 
corrections are enhanced by large logarithms which need to be resummed to all orders to obtain 
reliable predictions. For total broadening the leading logarithms have the form a" ln 2n (i?r). 
For thrust the resummation has been performed to next-to-next-to-next-to-leading logarithmic 
(N 3 LL) accuracy [3], while for broadening only NLL resummation is available [I]. For both 
observables, the next-to-next-to-leading order (NNLO) perturbative corrections are known 
[HIE]- For thrust, combining the fixed-order corrections with the resummed result, an accurate 
value for the strong coupling constant has been determined from a global fit to the available 
experimental data [7]. This yields one of the most precise determinations of a s currently 
available. Interestingly, the value comes out almost 4a below the world average [8]. It is 
important to validate this result using other event-shape variables. To this end, we extend in 
this paper the resummation for the broadening distributions to NNLL accuracy. 

So far, the only other event-shape variable for which the resummation has been carried 
out to higher accuracy is the heavy jet mass [9]. In the two-jet limit it is, however, closely 
related to thrust. In contrast, the jet broadening is not only literally orthogonal, but also 
its theoretical description in the two-jet region differs in important aspects. Both observables 
receive contributions from soft emissions and from radiation collinear to the jets, but the 
momenta scale differently in the two cases. For broadening, the transverse momentum of all 
emissions is of the same order as the broadening, \p^\ ~ \p^r\ ~ br- In contrast, the transverse 
momenta of soft particles contributing to thrust are power suppressed, {pjr/Pc) 2 ~ (1 — T) 
near the end-point T — 1. Because of this scaling, broadening is sensitive to soft recoil effects, 
while thrust is not. These recoil effects modify the higher-log resummation [I] and were missed 
in an early analysis |lUj . 
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For the above reason, also the description of the two event shapes in Soft-Collinear Effective 
Theory (SCET) [ITHT3~| differs. The effective theory relevant for thrust involves soft particles 
with power suppressed momenta, which are also called ultra-soft, while the components of the 
soft momenta are unsuppressed for broadening. The versions of SCET relevant in the two 
cases are sometimes called SCETi and SCET n . While the effective-theory analysis of SCETj 
observables is well understood, this was not the case for SCETn. One problem was that the 
loop integrals in the effective theory were not always well-defined in dimensional regularization, 
another that renormalization group (RG) invariance required the low-energy effective theory 
to depend on the hard momentum scale, which corresponds to the center-of-mass energy Q. 
This dependence arises naturally in SCETj from the power suppressed momentum of the soft 
modes, but is naively absent in SCETn. Both issues are now understood. The additional 
singularities can be regularized with an analytic regulator. The individual soft and collinear 
contributions then suffer from singularities in this regulator, which cancel once the contribu- 
tions are combined. However, the presence of the regulator breaks a rescaling symmetry which 
is not recovered when the regulator is removed. This effect is called the collinear factorization 
anomaly [14]. It induces Q-dependence in the low-energy theory. Consistency conditions can 
be used to show that this Q-dependence exponentiates [THEE]. Instead of deriving the ex- 
ponentiation from consistency conditions, one can also obtain it from an RG framework, the 
so-called rapidity RG [T6|fT7]. The associated framework is special in that one does not evolve 
from a higher to a lower virtuality, and that there is no physical coupling constant associated 
with the evolution. Nevertheless, using this formalism one arrives at the same exponentiated 
Q-dependence. The early papers on SCETn processes used versions of analytic regulariza- 
tion which break important properties of the effective theory such as gauge invariance and 
the eikonal structure of soft emissions. We have recently shown that only the phase-space 
integrals suffer from the additional divergences and that the regulator can be introduced in 
such a way that the factorization structure of the effective theory is maintained [IB]. The 
recent work [T5] and our paper provide the first examples of two-loop computations in SCET 
for anomalous observables. The simple form and the good properties of our regulator were 
crucial to make these computations feasible. 

Based on the new formalism, an all-order factorization theorem for the broadening dis- 
tributions was derived in [2U]- This allows for higher- log resummations, but in contrast to 
thrust the available perturbative input was only sufficient for NLL accuracy. In this paper, we 
derive the missing ingredients for NNLL resummation. These include the one-loop jet and soft 
functions as well as the two- loop anomaly exponent. The necessary computations are quite 
involved. In addition to polylogarithmic functions the results also contain elliptic integrals. 
As a check, we have independently performed the computations analytically and numerically. 
At small values of the broadening, our computation fully predicts the broadening distributions 
at 0(oig), and also all logarithmic terms at 0(a^). By comparing to numerical fixed-order 
results from event generators, we verify that we indeed obtain the correct behavior in the 
two-jet limit. Finally, we present numerical predictions for the resummed rates and discuss 
scale uncertainties and scheme dependence. 

Our paper is organized as follows. We start by discussing the factorization theorem for 
jet broadening and defining the relevant jet and soft functions. In Section [3] we perform the 
one-loop computations of the jet and soft functions. The main result of our paper is the 
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computation of the two-loop anomaly coefficient, which is described in Section HJ The two 
sections concerning the perturbative computations are necessarily technical. The reader only 
interested in the results can find these in Sections 13.31 and 14.31 We then derive the NNLL 
resummed result in Section [5] and present predictions for the broadening distributions at LEP 
energies. In Section we compare the perturbative expansion of the resummed result to 
numerical results of fixed-order event generators. The necessary expansion coefficients are 
listed in Appendix |A] The numerical evaluation of the jet and soft diagrams is discussed in 
Appendix [Bj and the relevant two-loop matrix elements are given in Appendix O Analytic 
results for the individual diagrams, both in momentum and in Fourier-Laplace space, can be 
found in Appendix ID1 

2 Factorization theorem for jet broadening 

Collider events with low values of the broadening consist of two nearly massless jets recoiling 
against each other, accompanied by soft radiation. The relevant factorization theorem has the 
form HSJEO] 

x J L (b L -b s L ,pi,n) J R (b R -b'to-p^n) S(b s L , b s R , -pi, -p^ n) . 

The hard function H(Q 2 ,fi) describes the production of the quark pair at large momentum 
transfer, before collinear and soft emissions. It is the square of the quark vector form factor 
H(Q 2 ,fi) = |Cy(— Q 2 — ie,fi)\ 2 , which is known to three-loop accuracy J2HI22]- What dis- 
tinguishes jet broadening from other dijet event shapes, such as thrust or jet masses, is that 
the transverse momentum of the soft emissions is of the same size as the one of the collinear 
partons. As a consequence, the jets recoil against the soft radiation, and the factorization 
theorem involves an integral over the transverse momentum. The broadening is the sum of 
the soft and collinear broadenings, and the transverse momentum of the collinear radiation is 
equal and opposite to the transverse momentum of the soft radiation. 

The broadening jet function for the left-moving collinear partons is defined as 

\ HU Mb L ,pi, AO = ^ (2n) d S(n ■ Px - Q) 5 d ~ 2 {pi - pi) 

X,reg. 

x ^ < X fe(°)l°> • ( 2 ) 

This definition involves the reference four-vectors n^ 1 = (1,Ht) and = (1, — n?) along the 
thrust direction. The field xi x ) is the collinear quark field of SCET. The collinear SCET 
Lagrangian is equivalent to the ordinary QCD Lagrangian and the field xi x ) can be identified 
with x( x ) — 4 W\x) ip(x), where ip(x) is the QCD quark field and W(x) a straight Wilson 
line along the n M direction from — oo to x (see e.g. [23] for more details). Also the soft SCET 
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Lagrangian is equivalent to the usual QCD Lagrangian and in all our computations, we will 
work with the standard QCD Lagrangian and Feynman rules. 

Beyond tree level, a naive definition of the jet and soft functions leads to integrals that 
are not fully regularized by dimensional regularization, because they suffer from rapidity di- 
vergences. The sum over intermediate states in the above definition contains a prescription 
to regularize these divergences. It consists in a modification of the phase-space integrals and 
amounts to the replacement [Tj 



Jd d k5(k 2 )6(k°) -> Jd d k(j±) 5(k 2 )6(k°), (3) 

where k + = n ■ k. As the rapidity divergences only occur in the phase-space integrals of the 
additional emissions, we find it convenient to use the standard measure for the tree-level phase 
space of the qq pair. The jet function Jr for the right-moving collinear partons is obtained 
by exchanging -h- in the definition (j2J), but keeping the same phase-space regulator 
(v + /k + ) a . Also the broadening soft function, 

S(b L , b R , P i, P i = ^ 5 d ~\pi L - pi) 5 d ~\pi R - pi) (4) 



X L ,X R ,reg. 
X 



i€X L j€X R 



is regularized in the same way. After the convolution of the jet and soft functions has been 
carried out, one can take the limit a — > 0. The individual functions suffer from divergences in 
this limit, which cancel in the convolution. 

The jet and soft functions, as well as the broadening cross section itself, are distribution- 
valued both in the broadening and in the transverse momentum variables. To work with 
regular functions and to turn the convolutions in ([TJ into a product, it is advantageous to 
perform a Laplace transform of the cross section 

db L e~^ / db R e-^-^— (5) 



dr L dr R J J db L db R 

and to work in impact parameter space rather than transverse momentum space. To this end, 
we Laplace and Fourier transform the jet and soft functions. For the left jet function, for 
example, we define 

/oo p Ad-2 _i_ 

db L e-^ J —J^e-**-* J L (b L ,pt,ri , (6) 

and analogously, for the right jet function and the soft function. One can further simplify the 
functions by averaging over the azimuthal angle, and define 

S{t l , t r , z l , z r , h) = J dtt%_ 2 J dQ%_ 2 S(T L , t r , xf , Xr, h) , (7) 
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where we have introduced the dimensionless variables Zl,r = 2\ x ~t r\/ t l,r- The normalization 
factor 

- ^d— 2 2 

= (27r) rf - 2 = (4vr) 1 - e r(l - e) (8) 

ensures that 5 = 1 at lowest order, for any value of the space-time dimension d = 4 — 2e. 
Introducing also 

^ L , fl (r, z,/j,)=J\f (27r) d - 2 - ( y ) Jl^t, ar 1 , //) , (9) 

the factorization theorem in Laplace-Fourier space takes the simple form 
1 d 2 a f°° f°° — — — 

— 3 j = H(Q 2 ,fl) dz L dz R J L (T L ,Z L ,fl)J R (T R ,Z R ,fl)S(T L ,T R ,Z L ,Z R ,fx). (10) 

cr dr L dr R J J 

As discussed above, the naive factorization theorems ([1]) and ffTOl) do not achieve complete 
factorization, since the additional regulator induces implicit Q- dependence. The main result 
of our previous paper [2D] was to derive the all-order form of this Q-dependence and to show 
that the product of jet and soft functions can be written as 

Jl(t l , z l , /i) Jr(tr, z r , h) S(t l , t r , z l , z r , h) 

= exp [-F b (t l , z L , n)\n (Q 2 fl) - F b (t r , z r , //) In (Q 2 f R )] W(t l ,t r , z l , z r , //) , (11) 

where the anomaly exponent F r (tl, zl, aO and the remainder function W(tl, tr, zl, Zr, /x) are 
independent of the large scale Q. The results (flOl) and (fTTT) form the basis of the all-order 
resummation of large logarithms. There are two sources of dependence on the large scale Q: 
the hard function H(Q 2 , fi) and the anomaly in (lllj) . The large logarithms in the hard function 
can be resummed by solving the RG equation for this function and evolving from a high scale 
\i>h ~ Q t° a lo w scale fi ~ 1/tl/r ~ \jh^j R . The logarithms associated with the anomaly 
are exponentiated in (fTTj) . To achieve NNLL accuracy, we need to evaluate the function 
W(tl,Tr, Zl, Zr, ji) to one-loop accuracy and the anomaly exponent F B (r,z,fi) to two-loop 
accuracy. The computation of the jet and soft functions to one-loop order will be discussed 
in the next section. The two-loop anomaly exponent can be extracted from a computation of 
the divergences of the two-loop soft function in the analytic regulator, as discussed in detail 
in Section H73l 



3 Jet and soft functions at one-loop order 
3.1 Soft function 

We first compute the one-loop soft function, which is simple since there is at most one soft 
particle in the final state at the given order. The one-loop virtual corrections are scaleless and 
vanish in dimensional regularization. The real emission diagrams are shown in Figure [TJ As 
the last two diagrams vanish and the first two diagrams give identical contributions, one is 
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Figure 1: Next-to-leading corrections to the soft function. The one-loop virtual diagrams are 
scaleless and vanish. 



left with a single diagram at this order. The result for the one-loop soft function has already 
been given in our previous work in a slightly different regularization prescription [20] . Setting 
— and v\ = v + Q, one can easily translate these expressions to the current scheme. In 
broadening-momentum space the result takes the form 



(12) 



a s C 



7T 



a 



Pl 



2-a 



6 b 



PL 



5(b R )5 d -\pi)-(L^R) 



with pi = \p~t\- After Laplace and Fourier transformation this yields 
S(t l ,t r ,z l ,z r ) = 1 



OI s Cf J 



i(I + 21n,f 
a \ e 



z 



8Li 2 (-^U41n 2 z^+ ; 



21 
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(L O R) 



where we introduced the notation z± = (-^/l + z\ ± l)/4 and t l = r L e yE . 



(13) 



3.2 Jet function 

The calculation of the jet function is considerably more complicated, since there are now up 
to two collinear partons in the final state, and the integrals involve a non-trivial angle in the 
transverse-momentum plane. At tree level, the left and right jet functions are given by delta 
functions, <JL R (b,p ± ) = 5(b — ^\p ± \), which after Laplace and Fourier transformation turn into 



Jl,r( t i z ) 



4 e T(2 - 2e) 

r 2 fi - 7) 



A-2e 



[1 + z 



2 \3/2-e 



(14) 



In the following, we will drop the subscript L, R of the tree-level jet function for convenience. 
At one-loop order, the left and right jet functions do not coincide anymore, since the regulator 
([3]) breaks the left-right symmetry, which is recovered only when the jet and soft functions are 
put together. The one-loop virtual corrections are again scaleless and vanish in dimensional 
regularization. Among the real emission diagrams, the last diagram in Figure [2] vanishes and 
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Figure 2: Next-to-leading corrections to the jet function. The one-loop virtual corrections are 
scaleless and vanish. 



the second and third diagrams give identical contributions. The jet functions can thus be 
written in the form 

Jl,r{^P L ) = J®^) + Jtn^p^ + 2jjM{b,p L ) + 0(a 2 ) , (15) 

where the term with the superscript (la) refers to the self-energy contribution in the first 
diagram in Figure [21 and the term (16) to the Wilson-line contribution in the second diagram. 
In the remainder of this section, we explain the analytic calculation of these diagrams. We 
discuss the left jet function in detail, but will only summarize the results for the right jet 
function. As a check of our calculations, we computed the jet functions numerically in Laplace- 
Fourier space. The details of our numerical approach can be found in Appendix IB. II 

Self-energy diagram 

In the original broadening-momentum space the self-energy contribution to the left jet function 
becomes 

J^\b^) = a s C F —(l-e)n 2€ j d d q5(q 2 )9(q°) J d d k (-±) 5(k 2 )9(k°) 

2k _ q 6(Q-q„- k_) o rf " V ~Q ± - k L ) S(b - \\q X \ - \\k^\) , (16) 



k_ 

x 



where q and k are the quark and gluon momenta, respectively, and fl 2 = /i 2 e 7B /(47r). Notice 
that we did not impose an additional regulator on the phase space of the quark, since we work 
in a scheme where the tree-level phase space remains untouched. We now use the first delta 
function in the second line to perform the /c_-integration, the second delta function for k 1 - and 
the on-shell conditions for the integrations over the plus components. We then arrive at 

o-l+2e rQ roc i>1 p 

J^ la \b,p ± ) = a s C F -^ r (l-€)^ J dq^ J dqq d ~ 3 J rfcos^sin d " 5 ^ J dQ d _ 3 

(Q-q-)v+ V Q-Q- 



^2 _j_ g2 _ 2pq cos 9 J q^p 2 + Q 2 q 2 — 2q_Qpq cos 9 

x 5(b - ^q - t^vV + Q 2 - 2pqcos9^j , (17) 

where 9 denotes the angle between q 1 - and p L , and we abbreviate q = Ig -1 ! and p = \p \. Next, 
we use the remaining delta function to perform the integration over cos#. It is furthermore 
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convenient to introduce dimensionless variables via g_ = Qrj, q = and y = p/2b. This gives 
r(ia) fh ±x _ / \ € C 1 - e ) n f \ % . 2s-i-e r (ia) , 



with the two-dimensional integral 

(19) 

Note that we factored out a distribution in (fl8i) . that is divergent in the collinear limit y — > 1. 

As the above integral is well-defined in the limit a — > 0, we may set the analytic regulator 
to zero. We then integrate over 77 and symmetrize the remaining ^-integral. Substituting 
6 = 1 + yy/u leads to the intermediate result 

(20) 

We next perform a partial integration and rewrite the resulting hypergeometric functions with 
their integral representations. After inverting the order of the integrations, the integral can 
be solved in closed form and we find 



1 - y 2 (1 - e)P(l - e) 



2(l- e )(l-y 2 ) 3 F 2 (l 1 i-e;|,l- e 



2' "'2 "'2'" 1-y 2 

1 3 _ 3 - 2 

2' '2 6; 2'" " , 1-y 



+ 2 / 2 (l-2e) 3 F 2 (-,l,--e;-,2-e;- T ^ ) ^ (21) 



In the collinear limit y — > 1, we derive 



#"(y i) - F Vrp-a") ° = 5 + (1 " ln2)£ + ° (f2) ' (22) 



As the prefactor of the integral in ()18|) diverges only for j/ -+ 1, the higher-order term is needed 
only in that limit. We may thus expand the exact expression for the integral as 

I ( L a \y) = \~ \V^7 2 ~ In (j^ + (1 - In 2)/(y) e + 0(e 2 ) , (23) 

where the precise form of the auxiliary function f(y) is not needed, but only its value at 
/(l) = 1. After taking the Laplace and Fourier transform we finally obtain 

x j-I + 2(l + 2 2 )lnz + -(2 + z 2 )ln(i±_^) -2^1 + ^ + 1 j, (24) 



with z + = (a/1 + z 2 + l)/4 and r = re 
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Wilson-line diagram 

The calculation of the second diagram in Figure [2] proceeds along the same lines, but this time 
it will be essential to keep the analytic regulator non-zero. We now start from 

jfW) = a S C F ^ e jfj d d q5(q 2 )6(q°) J d d k (j±Y 5{k 2 ) 9(k°) (25) 

?_(fc_ + 

k_ (2k ■ q) 



< q : {k - + q : ] 5( Q - q _- V - - k 1 ) S(b - \\q L \ - \\k- 



which results in 



a 



^(»,/) = %i (f-) (i - y'r 1 -' (26) 



with the two-dimensional integral 



LW ^ ^ y //l /; (£ -2^)2 + 477(1 -j/)(1 + j/-0 

(27) 

Because of the singularity of the integrand for rj — > 1 this integral would be ill-defined without 
the additional analytic regulator a. In physical terms this is the limit where the minus 
component of the gluon momentum, which we assumed to be of the order of the large scale 
Q, vanishes. We then expand in the analytic regulator using the relation 



[i-vr 1+a = - s(i-v) + 

a 



n 



+ 0{a) . {2i 



The subsequent calculation follows the same steps as before, but it is considerably more 
involved. In particular, the exact result in terms of hypergeometric functions now becomes 
rather lengthy. All we need here is the full expression in the limit y — > 1, as well as the leading 
term in the e-expansion for arbitrary y < 1. In the collinear limit, we now find 

rd^, ^ „ _oi-2e T 2 (2 - e) . _ , _ e . r(l-26)r(l + 2e) 



l£°>(y 1) ~ -2 1 -* ' + (1 - y) 



eT(3-2e) v * J r(l-e)r(l + e ; 

|I-l-ln(l-y)-31n2 + *(i + 6)-*(c)} 



(29) 



where \l/(x) = Y'(x)/Y(x) denotes the digamma function. In contrast to the self-energy dia- 
gram, our result for the jet function in f)26p now contains non-trivial distributions of the form 
(1 — y)~ 1_e , (1 — y)~ 1 ~ 2e and (1 — ?/) -1 ~ 2e ln(l — y), which shows that one cannot naively 
perform the e-expansion to extract the singular terms. There is no such complication for the 
terms with y < 1, where we obtain 



if>(v <i) = ^(2- ^) + 2(2 + yr^7) In ( i+^HZ ) 
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iy + y 2 



Hl-y)- 1 + 8 o y + y2 \n(i + y) + y /T=tf-l + 0{e). (30) 

We next combine the two representations in (I29p and (13"0"j) to account for the non-trivial struc- 
ture of the integral in distribution space. The subsequent Laplace and Fourier transformations 
can again be done analytically. Our final result for the Wilson-line diagram takes the form 

J L {t,z) = J {t,z)—^-{^t') [u+Qt ) (31) 
x j-i(i + 2mz+) + J^+ 2 -8Li 2 (-^) +8Li 2 (-Vl + 2 2 ) -41n 2 z + 

IT 2 > \ 

+ ln 2 (l + z 2 ) + 2z 2 ln(l + z 2 ) + 4(1 - z 2 ) ln(4^+) + 4Vl + z 2 - 8 In 2 - — \, 



IE 



where again z± = (a/1 + z 2 ± l)/4 and r = re 
Right jet function 

We now briefly summarize our results for the right jet function. First of all, we have seen 
above that the self-energy diagram does not generate a divergence in the additional regulator 
and, consequently, we obtain the same contribution as for the left jet function, 

J { n\T,z) = J { l a \T,z). (32) 
The Wilson-line diagram then follows from (126 j) and (j27p with the replacement 

"^-VWY^ (33) 



(2-0 2 b 2 J \(1-V)Q 

Note that the phase-space regulator now enters with the large scale, since k + ~ Q in the 
anti-collinear region. Without going into further details, the exact expression in broadening- 
momentum space becomes 

I™(v -+ 1) ~ -2 1 "* r2(2 " £) + (1 - v)" e 2 1 - r(1 - 2e)r(1 + 2£) 

^ u/-n;_ 2 er(3 _ 2e) + u yj ^ r(i-e)r(i + e) 

-- - 1 + ln(l - y) - In 2 - 2 7s - 2*(2e)l (34) 
a 2e J 

in the collinear limit. For y < 1 we get 

4% < 1) = ~ (2 - + l ~±f In (|^) + >/W - 1 + ^(e) • (35) 

After Laplace and Fourier transformation this translates into 
J r ( T ? z) = J (r, z) (ft t ) {v + /Q) 
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\-(- + 2\nz+) + - + 8Li 2 (-v / l + ^ 2 ) + ln 2 (l + z 2 ) + 2z 2 ln(l + 



- 2 ) 



4(1 - z 2 ) ln(4z+) + 4V1 + 2 2 - 8 In 2 



27T 2 



(36) 



3.3 Remainder function W 

We are now ready to combine our results for the jet and soft functions. According to (flip , 
the product of jet and soft functions determines the anomaly exponent Fb(t,z,/j,) and the 
remainder function W(tl, tr, Zl, zr, fJ>)- As the one-loop anomaly exponent has already been 
given in our previous work [20], we focus here on the extraction of the remainder function. 
The renormalized function W = Zy/W haie fulfils the renormalization group (RG) equation 

,) In (fi 2 f L f R ) - 4>y q (a s ) W(t l , t r , z l , z r , fi) , (37) 



W(t l ,tr,z l ,z r ,ii) 



dlo.fi 

where r cusp denotes the cusp anomalous dimension and 7 9 refers to the quark anomalous 
dimension as defined in [23]. The Z-factor Zw fulfils the same RG equation. To one-loop 
order, the solution takes the form 



'W 



47T 



-^-i(^ln( / , 2 f L f i? )-2 7( 



0{a s 



(38) 



where Tq = ACr and 7q = — 3CV are the leading coefficients in the perturbative expansions 

r cU sp = X^o r n(f^) n+1 and l q = J2n=o^fn (f^)™ +1 - To the g iven accuracy, it is furthermore 
convenient to rewrite the remainder function as a product of a left and a right function, 

W(t L , Tr, Z L , Zr, jl) = w(t L , Z L , fi) w(t R , Zr, fi) + O(o%) . (39) 

Combining our results for the jet and soft functions from the previous section, we find 

ln 2 (/i 2 r 2 ) + 31n(/i 2 r 2 ) - 8Li 2 (- ^) + 16Li 2 (- V / TT^ 2 ) 



W[T,Z,fi) = J [T, Z) < 1 + 



47T 



4 In 2 z + + 2 ln 2 (l + z 2 ) + 2(5 - 3z 2 ) \n(4z+) + 6Vl + z 2 



7T 



(2 - 3^)ln(l + zr) - 121n2 + 1 + — 



(40) 



The fact that our result for the remainder function W is left-right symmetric and that it fulfils 
the RG equation (f3Tj) provides a non-trivial check of our calculation. 

Later, we will need logarithmic moments of the remainder function. For the lowest mo- 
ments, we find 



dzw{r, z, fi) = 1 + 



a s CF 
4tt 



lnVf 2 ) + 31n(/r¥^ 



57T 2 

IT 



dz In z + w(t, z, fi) 



C^sCf 

Air 



6 + — -20C 3 + -7r 2 ln2- 121n 2 2- — ln 3 2 
2 3 3 



(41) 



These are used to compute the fixed-order expansion of the resummed results in Section [HI 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 

Figure 3: Diagrams that give non- vanishing contributions to the soft function at NNLO. In 
addition there are mirror-symmetrical graphs, which we take into account by multiplying each 
diagram with a symmetry factor Sj, where s a = Sb = Sf = s g = 2, s c = 1 and Sd = s e = Sh = 4, 
see text. 



4 Two-loop anomaly coefficient 
4.1 Setup of the calculation 

The anomaly coefficient Fb{t,z,h) can be extracted from the divergences in the analytic 
regulator of any of the soft and jet functions. Here, we will consider the two-loop soft function, 
since in this case there are again up to two partons in the final state and the integrals are 
similar to the ones that we encountered in the calculation of the one-loop jet function. As the 
divergences cancel in each hemisphere independently, it is sufficient to consider emissions into 
one of the hemispheres only. To be specific, we focus on the emissions into the left hemisphere. 

At two-loop order the purely virtual corrections are again scaleless and vanish. Among 
the mixed virtual-real and double real emissions, only the diagrams in Figure [3] give non- 
vanishing contributions. The same matrix elements also arise in other two-loop computations 
of soft functions [23H32], what makes our case different are the phase-space constraints and 
the necessity of working with an additional analytic regulator. In the following, we denote the 
individual contributions of these diagrams to the soft function by 

S?\b L ,b R ,pi,pi) } i = a,b,...,h (42) 

and similarly for the Laplace-Fourier transformed soft function S. 

The diagrams in Figure [3] have counterpart diagrams that follow by exchanging -h- 
as well as complex conjugation. It turns out that the matrix elements of diagrams (d), (e) 
and (g) are not symmetric under n M -H- n^. In an integral over a symmetric phase space, 
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these mirror graphs would nevertheless give the same result as the original diagrams. As 
our regulator fl3]) breaks this symmetry, it is however not guaranteed that the contributions 
from the mirror graphs and the original diagrams agree. We therefore explicitly calculated 
these mirror graphs and found that they do give the same results as far as the divergences in 
the phase-space regulator a are concerned. The contributions from these mirror graphs are 
therefore included in the symmetry factors Sj. 

Before going over to the more technical aspects of the NNLO calculation, we give the result 
for the one-particle cut in the last diagram of Figure |3j As there is only one soft particle in 
the final state, this diagram can easily be calculated along the lines of the NLO calculation. 
In broadening-momentum space, we find 

c {2h) ( , , j. u a 2 s , 2 , 2 e v% T 2 (l + e)r 3 (-e) 

Si '(b L ,b R ,p L ,p R ) = ^-C A C F (/i e 7B ) -± f(-27j cos ( 7re ) 

s(b L -f)6(b R )5 d - 2 (p^ } (43) 



x Pl 2 ~ 2e ~ a 



which after Laplace and Fourier transformation turns into 



S L {t l ,t r ,z l ,z r ) = y—J C A C F (fi t l ) 



a 



4 16 In 2+ 
^3 — 



f 48Li 2 (- Z -) + 32 In 2 z + + ™p) I + 288 S li2 (- ^ 



- 48Ll3 (" %) - " 192lnz+ Ll2 (- v) ~ nr lnz + ~ x C3 } (44) 

up to finite terms in the a-expansion that are not needed to extract the anomaly coefficient. 
The function Si^z) is a Nielsen generalized polylogarithm. Since we only consider emissions 
into the left hemisphere, the result only depends on zl and we therefore drop the superscript 
L on z± = (y/l + z\ ± l)/4 throughout this section for brevity. Next, we will discuss the 
calculation of the two-particle cut diagrams in detail. We will then explain how to extract the 
bare anomaly coefficient from the two-loop soft function, discuss its renormalization and give 
our final two-loop result. 



4.2 Two-particle cut diagrams 

For the two-particle cut diagrams, we only need to consider the case where both of the gluons 
are emitted into the left hemisphere. Similar to the one-loop jet function, the calculation 
is complicated by a non-trivial angle in the transverse-momentum plane. We organize our 
calculation along the same lines as in Section 13.21 It turns out, however, that one has to 
compute the two-loop soft integrals to one order higher in the e-expansion, which severely 
complicates the calculation. 

As an example of a two-particle cut diagram, we will consider diagram (b) in detail. The 
calculation of the remaining diagrams in Figure |3] proceeds similarly, and the results are 
summarized in Appendix [Dl As in the calculation of the one-loop jet function, we checked our 
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analytic results with a purely numerical approach. The details of the numerical calculation 
can be found in Appendix IB. 21 

We now turn to the evaluation of diagram (b). In broadening-momentum space we start 
from 



5(l 2 )6(f 



s? b \b L ,b R ,pi,Pi) = **c F (c F -^) fi 4e J d d k (j±ys(k*)9(k ) I d d l (?± 

(45) 

We now use the first delta function in the second line to perform the /-'--integration and the 
on-shell conditions for the integrations over the minus components. We then arrive at 



s~i ry — 1+4<e /-oo rK r 

st\b L ,b R :Pi,Pi) = o? 8 c F (c F - -f) A 4e y o dk k d - 5 y o dk + y o 

r 1 r f v 2 \ a 1 

x / dcos9 sin d_5 / dfi rf _. ' + 



' m ax 



i 

x r -^- - i * - ^/pi + fc2-2p L fccos6l) ^(Ph) W . (46) 



/ P/+ + (p| + A; 2 - 2p L & cos 9)k A 



with /™ ax = y/p 2 + k 2 - 2p L A; cos0, where # denotes the angle between /c- 1 and pf, and we 
write k = {k^l and pi = \p~t\- Next, we use the remaining delta function to perform the 
integration over cos#. Substituting l + = pk + , one may easily perform the /^-integration. We 
furthermore introduce dimensionless variables via k = and y = which results in 

S? b \b L ,b R ,pi,P^) 



(47) 



2tt 3 

with the two-dimensional integral 



A-y JO 1+P (2-0 

(48) 

Note that we again factored out a singular distribution in ( 14"Tj) . and that the remaining integral 
is finite in the limit a — > 0. The subsequent calculation proceeds along the same lines as in 
Section 13.21 In particular, setting a = we could again find an exact representation of the 
above integral in terms of hypergeometric functions. Again, we need the exact expression only 
in the limit y — > 1. We find 

if^iv -> i) - -(i - ,r VSgi:? {m - ,j - h2 - .(i + «) - . 

(49) 
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In addition, we need the first two coefficients in the e-expansion for arbitrary y < 1. To 
this end, we use the Mathematica package Hyp Exp [331 El], which allows one to expand 
hypergeometric functions around integer and half-integer parameters. The extension to half- 
integer parameters is currently restricted to special classes of hypergeometric functions. As 
some of the functions that enter our calculation do not belong to these classes (these are of 
the type 3 2 in the notation of [31]), we use an integral representation for these hypergeometric 
functions to perform the e-expansion manually. We could moreover make use of results from 
[35] for some special cases such as 3-F 2 (l, §, §; 2 + e, § + e; y 2 ). The expansion of the above 



2' 2 

integral for arbitrary y < 1 finally reads 

r (26)/ ^ ^ V 1 -y 2 ■ n 2n . JV ( l - \l l -y 2 \ t • fy /1 ~y 2 - 1 

r q ' (y < 1 = — arcsin y In 1 — y + < Li 2 iL ^=^= — Li 2 — , 

s K } 2y 2 1 y> \ 2 \i + ^i^i) 2 \i + ^/i^y2 



+ ^4C1 2 (7r — arcsin — Cl 2 — 2 arcsin J 



1 



- ln(l - y 2 ) ln(l + ^1^) + - ln 2 (l - y 2 ) je + 0(e 2 ) , (50) 

where Cl 2 (x) denotes the Clausen function 

Cl 2 (x) =-J dt\n (2sin(t/2)) = l -(hi 2 {e~ ix ) - Li 2 (e")) . (51) 

After Laplace and Fourier transformation this translates into 

c( 2b V \ ( as \ 2 n (n Ca \< 2-2\^i u + f L) 2a 
S L ( t l, tr, z l , z r ) = [^) C F \C F -—j {fi t l ) 

x || + ^± + (24Li 2 (- ^) + 16 In 2 z + + 3vr 2 ) ± + 64 /^l) - 80 S 1)2 (- ^) 

+ 56Li 3 (- Z -) + 961n^ + Li 2 (- Z -) + ^ In 3 2+ + 4vr 2 In (^) + ^Cs}, (52) 

which again holds up to finite terms in the a-expansion. The definition of the function h^{z) 
can be found in equation (1641) below. It involves elliptic integrals and cannot be expressed in 
terms of polylogarithmic functions. 



4.3 Extraction of the anomaly coefficient 

In our previous work [2D], we showed that, to all orders in perturbation theory, the logarithm 
of the soft function can be written in the form 

ln<S(r L , t r , z l , z R ) = 2F b (t l , z l ) \n(u + f L ) - 2F B (r R , z R ) ln(v + f R ) + k% (r L , r R , z L , z R ) , (53) 

where we suppressed the divergent terms in the analytic regulator that cancel in the product 
of the jet and soft functions. The function k$ contributes to the remainder function W and is 
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therefore irrelevant for our purposes. In our two-loop calculation, we restricted our attention 
to emissions into the left hemisphere and we focused on the divergent terms in the analytic 
regulator. We will now show that this is sufficient to extract the anomaly coefficient Fr{t, z). 

Making the divergences in the analytic regulator explicit, the logarithm of the soft function 
can be rewritten in the form 

oo - 

In S(r L , t r , z L , z R ) = ("+TL) na (- f n (r L , z L ) + r n ,(r L , z L j) +..., (54) 

ra=l 

where the ellipsis denote terms that involve the variables of the right hemisphere. According 
to our phase-space regulator (jSJ), an n-particle cut diagram involves a power (v + ) na . The 
functions /„ and r n are therefore of G(ot™). Expanding the above decomposition to two-loop 
order and keeping only terms that are relevant for the extraction of the anomaly coefficient, 
we arrive at 

h\S(r L , t r , z l , z r ) = (7i(r L , z l ) + f 2 (r L , z L )^j - + (/i(t l , z l ) + 2f 2 (r L , z L ) \ hx(v + f L ) + ... . 

(55) 

Notice that the coefficients of the divergence and the logarithm do not agree beyond one-loop 
order. This is different from the more common situation in dimensional regularization, where 
the renormalization scale is tied to the coupling constant in the form a s fi 2e . In contrast, the 
two- loop diagrams in our calculation do not have a homogeneous scaling in u + , since they 
involve both one-particle and two-particle cuts. 

From fl53|) we see that the anomaly coefficient is determined by the logarithm rather than 
the divergence. Up to two-loop order it is given by the combination 

F B (r,z) = ^f 1 (T,z) + f 2 (r,z). (56) 

Taking the exponent of (I54|) and expanding to second order in a s , we may then extract the bare 
anomaly coefficient from the single logarithmic terms in v + . As this also generates a crossed 
term between the one- loop terms in j\ and r±, we first need to determine their higher-order 
terms in e. We find 

h(r, z) = ^(AT [l + 21n^ + + (2Li 2 (- + 21nV + £)e + (W- ^) 

- 4 S lj2 (- ^) + ~ In 3 z + + 4 In z + Li 2 (- ^) + y In z + + ^Ca) e 2 + G(e 3 ) } 

(57) 

and 

n(r, z) = ^(A 2 ) £ {" | + 8Li 2 (- Z -) + 41n 2 z + + ^ + (sLi 3 (- *-) (58) 
- 32 S ll2 (- Z ~)+f In 3 z + + 24 In z + ^(~ Z -) + ^- In z + + |cs) e + Q{e 2 ) } . 
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We next renormalize the coupling constant in the MS scheme, and recall that the renormalized 
anomaly coefficient fulfils the RG equation [20] 

d 



dlia.fi 



F b (t,z,ii) = r cuS p(a s 



(59) 



The anomaly coefficient renormalizes additively, Fg = F^ am + Zp B , and the Z-factor fulfils the 
same RG equation. Up to two-loop order, the solution takes the form 



J F B 









rrf 




4tt 2e " 


V47T/ 


) 


4e 


4e 2 



0(al 



The relevant anomalous dimensions are 

lt = 4C F , = ["--■-)C A -^T F n f , 



Tf /67 n 2 \ „ 20 „ 



1 



9 



9 



(60) 



(61) 



The renormalized anomaly coefficient finally takes the form 



a 



F b (t, z,fi) = — T { ]n(z + fiT) 



Air 



(3 \n 2 ^f) 



rf 

r, 



, 2(3 \nz + ) \n(fif) + d B (z) 



(62) 

with z + = (vl + z 2 + l)/4. The nonlogarithmic piece of the anomaly exponent is encoded in 
the function d B {z), for which we find 



d B (z) = C A 



1 + z 2 



h{z) + 



67 + 2z" 



h 2 (z) -8/13(2) + 32S 1>2 (- ^) - 8Li 3 (- ^) 



1 — U> 



9 v ' 9 ' v 2+ 

+ 8 S h2 (-w) - 24Li 3 (-u;) - 24S 1)2 (1 - w) + 8Li 3 (l - w) + 24S lj2 

-8Li 3 (^-^) -8(31nz++41n2) Li 2 (-^) +8 In ((1 + wW 3 ) Li 2 (-w] 



- 8 In 2 Li 2 



11 



l — io 



in 1 n 

4 In - ln 2 z+ + 121n 2 u? ln(4z+) hi 3 (2z+) 

3 



+ — ln z ^+ + 161n2 In — lnz+ + (24hT 2 
3 

+ it 2 In 2 + 
2(1 + * 2 



67 
~9~ 



In z + + 4 In 2 2 In 



w 



^ - 18(3 - !*> 
27 S3 9 



2w(32-z 2 ) 1 /1 + «A w(65 + 2z 2 ) 



9 



In 



w 



9 
4 



9 



, . , 2 13 + 22 2 , . . 4. 2 20, 4 2 82 

-/i 2 (^) In z+ \nz + + -z 2 

9 v ; 3 + 9 + 9 27 



4w(5-2 2 ) 1 /1 + un 2w(ll + 2z 2 ) 
H ^ -In ' 



9 



w 



9 



(63) 



where w = a/1 + z 2 and z± = (w ± l)/4. The result involves the three functions 



arcsin i 



Vl + t 2 z 2 
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/ , arcsint , r-— . , 

h 2 (z) = I dt -^===^lT¥^ : (64) 



o 



h 3 (z) = f dt ^2^= In (1 + VTT¥z^) 
Jo Vl — t 2 



which cannot be expressed in terms of polylogarithmic functions. They are related by the 
differential equations 



n 2 



z h' 2 (z) = h 2 (z) - hi(z) , z h' 3 (z) = — - hi(z) , (65) 

o 

and can be written in terms of elliptic integrals, 

hl (z) = ^F^,-z 2 )- £ d9F(9,-z 2 ), 



7T 



h 2 {z) = -E(l,-z z )- d9E(9,-z 2 ), (66) 



o 



where 



r e 1 f 6 i 

F(6,x)= / dt = , E(9,x) = / dt VI- xsin 2 t 

Jo y 1 — x sin 2 1 Jo 



(67) 



denote the incomplete elliptic integrals of the first and second kind, respectively. For 9 = tt/2, 
the functions are also called complete elliptic integrals and denoted by K(x) = F(^,x) and 
E(x) = E(^,x). The integrals involving the elliptic integrals in ( |66|) can be expressed in terms 
of the Kampe de Feriet hypergeometric function [56] . 

For the numerical implementation of the resummation, the expansion of the elliptic func- 
tions around z = oo is useful. We find 

, . , 2G 1 1-2G 4 

hl{z) = --7 2+ ^z^ + 9z^ + --- 

j z \ G 1 1 -2G 4 

h ^ ) = Z + --^ 2+ ^ + 4^ + -- (68) 

, , , vr 2 , 7 _ 7T 2 , 2G 1 1 - 2G 1 

h 3 (z) = — In z H C3 hi 2 H H — H + . . . . 

A ; 8 16 s 8 z 2^ 2 12z 3 9z 4 

The expansions involve Catalan's constant G ~ 0.915966. 

Finally, for the fixed-order expansion of our result, we will need the integral of the product 
of the tree-level jet function with the anomaly exponent. We find 

00 (0) j UG 34 81n 2 2 201n2 
dzj (z) d B {z) = T F n F I — + — — + — — 

„ f43G 17C 3 104 16, , 22 In 2 2 64, 2 , 1 
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5 Resummation 



Our result for the cross section in Laplace space has the form 



(fa 



POO POO 

h(q 2 ,v) dz L ^(g 2 fl)- Ffl(Li ^^(g 2 ^)- FB(ta ^ 

Jo Jo 



(To dr L dr R J J (70) 

x W(L l ,L r ,Zl,z r ,h) . 

For later convenience, and with a slight abuse of notation, we have written the anomaly 
exponent Fb as well as the remainder W as functions of the logarithms 

L L = \n(pf L ) and L R = \n(/jf R ) (71) 

instead of tl and r R . To get the resummed result, we need to solve the RG equation for the 
hard function and perform the Mellin inversion of the above result. 

To perform the inversion, we first rewrite the anomaly function in the form 

F B (L, z, /i) = ^r F [L + \nz + + f B (L, z, fj)] , (72) 



where z + = (vl + z 2 + l)/4 and the function 



f B (L,z,n) = 



/3 L 2 + ( J ± + 2f3 ]nz+)L + d*(z) 



0{al) (73) 



V = Vf = ^r^ln — . (75) 



collects the higher-order contributions to the exponent. The anomaly can then be written in 
the form 

^ Q2f2 yF B (L, z ,») = ^ f Y^ z+ Y^ e -2LF B[ L, Z )^ FfB{ L, Z ) (?4) 

with 

—Tq In — 

Distinguishing rj and rjp in (174|) will allow us to rewrite powers of the logarithm L as derivatives 
with respect to 77. Since the coupling constant is multiplied by a large logarithm for fi ~ 1/r, 
the quantities 77 and 77^ count as 0(1) and must be kept in the exponent. The remaining 
exponential on the right-hand side of (!74|) can be expanded in a s , since it does not contain 
large logarithms. After this expansion, the r-dependence of the cross section becomes very 
simple: at any fixed order in a s , it is given by a polynomial in the logarithm L multiplying the 
factor (fif)' 71 . The Mellin inversion can then be evaluated in closed form, using the formula 

5s Lj^^ = { - a " r l © m ■ (76) 

Performing the inversion for the cross section, we obtain the following result for the double 
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differential cross section 
1 d 2 



■'a 



W(-d VL ,-d VR ,z L ,z R ,n) 



x exp [2d VL F B (-d VL ,z L , /i) - rj F f B (-d VL ,z L , n) + (L -B- R)] 



1 fb T .\ VL e~^ L 1 



x 



b R \ m e -^ EriR 



r(r/i) b R \jj, ) T(rj 



R 



Vl=Vr=Vf 



(77) 

Since the exponential in the second line can be expanded in a s only a finite number of deriva- 
tives need to be computed, when evaluating the cross section at any given accuracy. At NNLL 
accuracy only the first and second derivative will be needed. 

To NNLL accuracy, the general expression can be simplified further by noting that the 
remainder function can be written as a product of a left and a right function 



sec 



W(L L , L R , z L , z r , ji) = w(L L , z L , /i) w(L R , z R , fi) + 0(a 2 s ) , 
. To obtain compact expressions one can then define the integral 

poo 

T(-L,fi)= / dz L (z+) ^ exp[-2LF B (L,z,n) - rj F f B (L, z, /jl)]w(L, z, n) 
Jo 



(78) 



(79) 



in terms of which the double differential cross section reads 

ji ) T(r] L ) b R \/i 



- J^r- = h(q 2 , *0 Ad VL ^) mm y) r ' 

a db L db R b L 



b R \ VR e^ Ef]R 



VL=VR=VF 



(80) 



For the total broadening 6 T = bi + b R and the wide broadening b w = m&x(br,,b R ), the 
expression simplifies further and we obtain 



1 da 



o"o dbx 
1 da 



H{Q 2 ^)X 2 {d r)T ^) 



bx 



'IT 



a db 



H(Q 2 ,n)I(d VL ,fi)I(d VR ,n) 



r) T =2r)F 

(Vl + Vr) e~~ lEivL+TlR) 1 

r(i + 77 L )r(i + 7te) b^ 



b w 



VL+VR 



Vl=Vr=Vf 



With the result for the cross section at hand, all that is left is to resum the logarithms 
in the hard function. This step is by now standard, since the same hard function appears in 
many processes, in particular in all dijet e + e~ event shapes and in Drell-Yan production at 
hadron colliders. The hard function is given by the square of the on-shell quark form factor, 
H(Q 2 ,fi) = \Cy{— Q 2 — ie, /i)| 2 , which is known to three loops j5T][22]. I n the following, we 
will suppress the ie-prescription attached to Q 2 . At one- loop order, the form factor is given 
by 



C v (-Q 2 ,fi h ) = l 



C F a s (/J, h ) 
4it 



7T 



-1/ + 3L - 8 + — 
o 



(82) 
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Figure 4: NLL and NNLL predictions for the total and wide broadening at Q — Mz- As a 
reference, we also plot the experimental results of the ALEPH collaboration. 



where L = \n(—Q 2 /fi\). The solution of its RG equation reads 

_Q2 \ - Q r(Mh>M) 



C v (-Q = exp [2S(fi h ,(i) - 2a 79 (/Xh, A*)] \ 



C v (-Q 2 ,fi h ), (83) 



where the integrals over the anomalous dimensions are defined as 



^ln — rf a r (/io,/i) 

(10 ^ 



7^ r L P ("s(/i)) 



i4) 



and similarly for the function a 7 g. Explicit results for these functions, together with the 
necessary anomalous dimensions can be found, for example, in Appendix B of [T4"] . 

In Figure HI we show the result obtained from a numerical evaluation of the differential 
broadening cross sections at NLL and NNLL. We normalize to the total hadronic cross section, 
whose perturbative expansion to NLO reads 



a = a 1 



Air 



55) 



The NLL and NNLL distributions are normalized to the fixed-order total cross sections at 
LO and NLO, respectively. We use dimensionless variables B T = bx/Q and B w = bw/Q, set 
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Q = Mz and use the PDG value a s (Mz) = 0.1184 for the strong coupling constant [8]. The 
upper two plots show the total broadening distribution, the lower ones the wide broadening. 
The bands in the plots are obtained by varying the scales by a factor two around their default 
values fi = bx (or fi = by/ for wide broadening) and fi^ = Q. The resulting variations 
are added in quadrature. The dominant uncertainty comes from varying the low scale /i. 
At NNLL (NLL), the hard scale variation is around 1% (15%). To use our results for a 
determination of a s , one will have to also include the matching to fixed-order results to 
obtain reliable predictions at larger values of Bt or By/- Furthermore, one will need to get 
a handle on hadronization effects, which modify the distributions in the peak region and 
below. Without accounting for these two effects, a comparison to experimental data is not 
completely meaningful, nevertheless we include the ALEPH data [3B] in the plots as a reference. 
Experimental results are also available from DELPHI [39J, L3 jlQ] and OPAL [4"T] . 

The plots on the left and right side in Figure H] show two different ways of setting the scale 
fi. We can either directly compute the differential cross section using (IHIl . or we can obtain 
it as a derivative of the integrated cross section 

R{B T )= f T dB T -^-. (86) 

Jo a dt >T 

For a fixed value of the scale fi it is trivial to integrate f l8"Tj) to obtain a resummed result 
for R(Bt). The only reason that the numerical results on the left and right side in Figure H] 
differ is the choice of the scale: choosing fi oc b-r = Q Bt in R(Bt) does not commute with 
taking the derivative with respect to Bt- The differences between the two prescriptions are 
not negligible. In particular, we observe that the scale uncertainties come out much smaller 
when working with R(Bt)- However, in both cases, we observe that the NNLL bands nicely 
overlap with the NLL result. 

In the traditional resummation literature pflfTO] , the differential distribution is obtained as 
a derivative of R(Bt)- An argument in favor of using this prescription is that the distribution 
is then properly normalized once the matching is included, since the integral over it produces 
R{Bt) and the resummation effects become small at large values of Bt- In fact, in addition to 
performing the resummation for R(Bt), the resummed logarithms are usually modified such 
that they vanish at the end-point i?™ ax , for example by replacing [I] 

w T ^w T ~i^ + 1 - (87) 

Once this is done i?(i?™ ax ) = 1. On the other hand, from an effective field theory standpoint, 
it is more natural to compute the differential spectrum directly (if this is the quantity of 
interest) and choose the renormalization scale according to the physical scales in the problem. 
Determining the spectrum by differentiating R(Bt), one effectively takes the difference of two 
large quantities, which are evaluated with slightly different scale choices. It has been argued 
that this leaves spurious contributions in the difference and should be avoided [7J. 
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Figure 5: Numerical comparison to results obtained using the Event2 generator. The upper 
left panel shows the coefficient B, as computed by Event2, and the lower panel the difference 
AB to our result. The difference must go to zero for large negative values of the logarithm 
L = In Bt- The right panel shows the contribution of individual color structures in AB. 



6 Fixed-order expansion 

It is instructive to perform the fixed-order expansion of the resummed result. To this end, we 
write the expansion in the form 
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The coefficients A and B can be computed numerically using the Event2 generator 
The NNLO coefficient C has been obtained in [5j[6]. Note that the expansion coefficients are 
normalized to the tree-level rate. In addition to the differential spectrum, we also consider 
the integrated rate R{B) for B = Bt or B = Byy, which was defined in ( 186]) . We write the 
expansion of the logarithmic part of R{B) in the form 
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The quantity Rl(B) is the leading term in the expansion of R(B) around small B. In the 
following, we will give explicit results for the expansion coefficients and will compare our 
result for Rl(B) with numerical evaluations of R{B) at small values of the broadening. A 
NLL computation determines the one-loop coefficients Gu and Gu as well as the two-loop 
coefficients G23 and G22- Our computation gives, for the first time, also the coefficient G 21 . 
For the total broadening, we obtain 
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In addition to ^-values, the result also contains Catalan's constant G ~ 0.915966. To obtain 
the coefficient, one needs to evaluate moments of the remainder function w(r,z,fi) and the 
two- loop anomaly exponent Fb{t, z, /x). These were given in equations (|4ip and (16"9"|) . With 
the constant G 2 i in place, the two- loop coefficient B(Bt) is known analytically in the limit of 
small Bt and can be compared to the numerical results of Event2. This is done in Figure 
[5j where we show the coefficient as a function of L = In Bt- At small values of Bt, the 
difference between the logarithmic terms and the full result must go to zero. We find that 
this is indeed the case within the numerical uncertainties, see Figure Instead of comparing 
the analytical result to Event2, we can also use the program to extract the coefficient G 2 i 
numerically. The values obtained in this way have an accuracy of about two per cent and 
agree with our result fl90l) within their numerical uncertainty, which supports the correctness 
of our computation and the validity of the resummation formula. The same is true for wide 
broadening; the relevant expansion coefficients are given in Appendix [A] Our values for G 2 i 
also agree well with the numerical results obtained earlier in [4]. 

Using our analytical results, we can extract the constant C2, which is the last unknown 
two-loop ingredient. To do so, we compute R(B) in the form 



R(B) = 1 



B 



a dB 



(91) 



The integral on the right-hand side of (|9~T|) can be evaluated numerically using Event2, taking 
into account the normalization to the total cross section a to NLO, which was given in 
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Figure 6: Determination of the two-loop constant C 2 for total broadening. The first plot is 
the full coefficient, the remaining ones show the individual color structures. 



We then consider the difference between the full result R(B) and its logarithmic part Rl(B). 
Since everything except the constant C 2 is known, the constant immediately follows from the 
requirement that the difference must vanish in the limit B — > 0. Figure E] shows the difference 
as a function of the logarithm as well as the value of C 2 we extract from it. For the total 
broadening, we obtain 

C 2 = (105.9 ± 0.5) Cf c a + (17.0 ± 0.4) C 2 + (-20.8 ± 0.3)c F T F n f = 102.3 ± 0.9 , (92) 

and the result for the wide broadening is 

C 2 = (128.5 ± 0.6) CfCa + (-0.41 ± 0.27) C 2 + (-15.0 ± 0.2) CpTFnf = 113.2 ± 0.4 . (93) 

The number for the total contribution has been obtained by fitting the sum of all color struc- 
tures. Within uncertainties, it agrees with the result obtained from adding the fit results of 
the individual color structures. Figure E] makes it clear that it is somewhat delicate to extract 
these numbers. On one hand, one wants to make the broadening as small as possible to avoid 
power suppressed terms which contribute to R(B) — Rl{B). On the other hand, one cannot 
make it too small because the numerics become unstable at very low B. To be able to reach 
very low values of B, we have run Event2 at quadruple precision and have generated 10 8 
events with a very low cutoff of 10~ 16 . To obtain the above values and uncertainties, we have 
performed \ 2 fits using all fit intervals in the range L G ( — 11.8, —4.6) containing at least six 
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Figure 7: Comparison with EERAD3 [5] for total broadening for the individual color struc- 
tures. The EERAD3 code uses an internal infrared cutoff yo f° r numerical stability. The three 
histograms correspond to yo = 10~ 5 (green), 10~ 6 (red) and 10~ 7 (blue). Our result is the 
black line. Larger values of the cutoff yo lead to deviations at small Ib.(Bt)- 



fit points. We then select the 12 fits with the lowest % 2 -value and use the spread among these 
as an estimate of the systematic uncertainty due to the choice of the fit range. The statistical 
and systematic uncertainties are of similar size and have been added in quadrature for the 
error estimates in fl92|) and (l93]h 

It is also interesting to compare our results to the numerical NNLO predictions obtained 
from EERAD3 [5]0 These provide the coefficient C(Bt) in flHHl) . which is plotted in Figure [71 
The comparison serves two purposes: first of all, it allows us to check that the resummation 
formula, by verifying that we correctly obtain the numerically dominant 0{a z s ) terms at small 

1 We thank Thomas Gehrmann for providing us with the results. 
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Figure 8: Same as Figure [7J except that we plot the difference between EERAD3 and our 
prediction for the singular terms. The difference must become constant at small broadening. 
Up to cutoff effects, this is indeed the case. 



values of the broadening. Secondly, it tests the numerical stability of EERAD3 at small values of 
the broadening and the sensitivity to the internal cutoff. To avoid numerical instabilities, the 
generator imposes a cut on the phase-space variables [13] . The observables are independent of 
the value of the cut parameter yo, if it is chosen small enough. By comparing to the numerical 
results at very small broadening, we can check the residual sensitivity to the choice of yo. We 
show three different histograms in Figure [7J corresponding to three different choices of y$. For 
the larger values of y one observes dramatic deviations at small B?, while we observe nice 
agreement for the smallest choice yo = 10~ 7 of the cutoff parameter. 

Our NNLL computation does not fully determine the logarithmic part at NNLO. It allows 
us to get the coefficients G34, G33 and G32, which are given in Appendix |A] but the coefficient 
G31 is beyond the accuracy of our computation. As a consequence, the difference AC(Bt) 
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between our result and the numerical evaluation using EERAD3 does not vanish, but becomes 
constant at small broadening. We plot the difference AC(Bt) in Figure [HJ Within numerical 
uncertainties the expected behavior is indeed observed, except at very low Bt, where the cutoff 
effects are non-negligible, as can be seen by comparing the results obtained with different 
cutoffs. We have performed the same numerical comparison for wide broadening and find 
consistent results also in this case. 

7 Conclusion 

In this paper, we have extended the resummation for the jet broadening e + e~ event-shape 
variables to NNLL accuracy. Our analysis is based on the all-order factorization theorem we 
obtained earlier [20J. The factorization for broadening suffers from a collinear anomaly, and 
to obtain NNLL accuracy, we have computed the anomaly exponent to two-loop order and 
the jet and soft functions to one- loop accuracy. 

The presence of the anomaly makes the computation subtle, since one needs to introduce 
an additional regulator beyond dimensional regularization to make the individual jet and soft 
functions well-defined. We recently conceived a convenient way of introducing the additional 
regularization, which is based on a modification of phase space, and respects gauge invariance 
and the factorization structure of the effective theory [18] . An alternative scheme was proposed 
in [T(3|ITT]. Together with the recent work [TJJ], our paper provides the first example of a full 
two-loop computation in SCET for an observable which is sensitive to soft recoil and the 
results demonstrate that our regularization method works in practice. 

The anomaly coefficient was extracted from the two-loop broadening soft function com- 
puted at fixed recoil momentum against the collinear radiation. The associated constraint 
on the two-emission phase space is quite complicated. Despite the complexity of the relevant 
phase-space integrals, we have managed to perform the computation analytically. In addition 
to polylogarithmic functions, we find that the result also involves elliptic integrals. As a con- 
sequence, the expansion coefficients for jet broadening do not only involve (^-values, but also 
/3-values (where /3 is the Dirichlet ^-function) such as Catalan's constant. Other examples of 
perturbative computations involving elliptic integrals can be found for instance in jmH3]. 

In addition to numerically verifying our analytical results for the diagrams, we have also 
compared the expansion of our resummed result in a s to the predictions from fixed-order event 
generators. Running Event2 with high statistics and a low cutoff, we have checked that we 
reproduce the numerical result for the wide and total broadening distributions at 0(a 2 s ) at 
small values of these parameters. Comparing to EERAD3, we have performed the same check 
at 0(af). Also in this case, we observe the expected behavior. 

Determinations of the strong coupling constant from LEP event-shape data typically use 
about half a dozen different event shapes, among them thrust, heavy jet mass, total and wide 
broadening [38THT| H6 | H7]. Statistically, there is little gain in using several observables since 
all of them are based on the same collider events. However, the consistency among the values 
extracted from the different observables provides an important cross check on systematic effects 
and on the reliability of the theoretical uncertainty estimate. With our results for broadening, 
higher-log resummation is now available for all four of the event shapes mentioned above. It 
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will be interesting to perform an extraction of a s based on these results and to check whether 
the low value of a s obtained from a global fit to thrust data based on N 3 LL resummation is 
confirmed [7j. 

Acknowledgements: We are grateful to Xavier Garcia i Tormo, Thomas Gehrmann, Tobias Hu- 
ber, Christian Lorentzen, Matthias Neubert and Gavin Salam for useful discussions. The work 
of T.B. is supported by the Swiss National Science Foundation (SNF) under grant 200020- 
140978 and the Innovations- und Kooperationsprojekt C-13 of the Schweizerische Univer- 
sitatskonferenz (SUK/CRUS). 



A Expansion coefficients 

Here we give the NNLL expansion coefficients defined in (1891) up to three-loop order. For total 
broadening, the coefficients read 
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The coefficient G21 was given in f l90|) . The three-loop coefficients are 

G34 = —^-CACrnfTp —C a Cf — —Cpn^Tp , 

+ ^ + 256^2 + W j c>iTf + ^ _ 128 b s 2 + 64 l2 fc 2 y r _ (95) 

G 32 = -41.46 C A C F n f T F - 200.62 C A C| + 78.75 C\C F + 97.09 C>/T F 
+ 9.103 C F n)T 2 F - 166.92 . 

We only list the numerical value for G 32 , but by computing higher logarithmic moments of 
the remainder and the anomaly functions, it could be obtained analytically from our results. 
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For wide broadening, the leading one-loop coefficients as well as the leading-logarithmic 
higher-order coefficients agree with the total broadening. For the remaining coefficients, we 
list the difference to the total broadening coefficients, for which we obtain 

AG 2 2 = , 
AG 21 = (32C 3 - 4vr 2 ) C 
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AG S3 = l -^C A Cl - ^Cl„,T F - 128( 3 Cj , (96) 
AGs2 = + 352C3 - ^) C A C- + (-128& + t> /3 > 

.^ + 288f3 + ^W2)cJ. 



B Numerical evaluation of the jet and soft functions 

As a check of our results, we have independently evaluated the diagrams numerically. The 
numerical evaluation is done directly in Laplace-Fourier space. The relevant phase-space 
integrals suffer from soft and collinear divergences, which we extract using a combination of 
subtraction and sector decomposition methods. 



B.l Jet function 

We start with the self-energy diagram and compute its Laplace and Fourier transform which 
yields 

~ 2 1+2e 1 

J { L a \r,x ± ^) = a s C F -^- e {l - e)/i 2e ^-p-^ J(r,x ± ,/i) , (97) 
where the relevant phase-space integral is given by 

J(r, x ± , n) = J d d q S(q 2 )9(q°) J d d k 5{k 2 )9{k°)5{Q - q_ - k. 
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(98) 

and we use the notation k = \k±\ = y— k\ and q = {q 1 ]- This is the same as expression 
(1161) . except that the delta functions have turned into exponentials because of the Laplace and 
Fourier transforms. 

The jet function only depends on the magnitude of x±, but not its direction. We can 
therefore average over the direction of x± relative to the total momentum k± + q± and replace 

-i(k ± +q ± )-Xj_ . d-3 / , , • d-4 1 „-i|fc ± +g ± \x cos <j> 
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inside the integral, with x = \x |. We then use light-cone integration variables and integrate 
over the k + , q + and g_ components, which leads to 

q2 pQ poo poo pit pit 

I(t,x ± ,h) = / dk_ / dkk d ~ 3 dqq d ~ 3 d9sm d ~ A 9 (fy sin*- 4 </> 

4 Jo Jo Jo Jo Jo 

■7j(k+q)+ixy/k 2 +q 2 +2cos9kqcos<f> (100) 
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with k + = k 2 /k_ and q + = q 2 /q~ and q^ = Q — k^ due to the three delta functions in fl98l) . 
We now eliminate x in favor of the dimensionless variable z = — and introduce new 

r 

integration variables u, v and r\ (u = 1 — u and v — 1 — v): 

k_ = Qu, k = —, q=—. (101) 

ru ru 

After the variable change, the integration over 77 can be performed, which yields the result 
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Without the term on the second line, we can perform the integral over v and 9 analytically 
in d = 4 — 2e dimensions 
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The divergence arises from the collinear singularity at the point v — ~ and 9 = 0. This is the 
only singularity of the original integral, and we can isolate it by rewriting it in the form 



I(T,x\ri = J Q dv jT d9 sin d -H v2 - V _ 2 ~_ V 2v _ - { [f(v, 9) - /(I, 0)] + /(I, 0) } . (104) 

The subtracted part in square brackets is finite for d = 4 and can be evaluated numerically. 
In the remainder, the w and 9 integrations can be evaluated in d dimensions using fllQ3p . Once 
this is done, the integral can be expanded around d = 4 on the level of the integrand and the 
remaining integration can be performed numerically. 

The jet-function diagrams involving Wilson-line emissions suffer from soft divergences in 
addition to the collinear singularity present in the self-energy diagram. In the parameterization 
(jlOip . the soft singularities appear when u and v go to zero and are regularized analytically. 
To isolate the two types of divergences, we split the v integration as 

1 pV cu t pl 

dv= dv+ dv. (105) 
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Figure 9: Sector decomposition for the numerical evaluation of the two-loop diagrams. Along 
the boundary of the integration region (red square) soft singularities occur, while the collinear 
singularities are located at v — | (blue dashed line). We cut the integration into eight patches 
along the dotted lines. 



As long as < v cut < 1/2, the first integral contains only soft divergences. The second part 
only suffers from the collinear divergence and can be evaluated numerically using the same 
technique as in the case of the self-energy integral. To compute the soft part, we first rescale 
v so that it runs again from ... 1 and then use sector decomposition [HH50] and split the 
integration over u and v into the sectors u < v and u > v. After rescaling u — > v u in the first 
sector and v — > u v in the second, the divergences are disentangled and can be extracted using 
relations such as 
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B.2 Two-loop soft function 

The one-particle cut diagram (h) in Figure [3] can easily be evaluated analytically directly 
in Fourier-Laplace space. For our numerical check, we therefore focus on the two-particle 
contributions. As in the jet function case, the corresponding integrals suffer from soft and 
collinear divergences. 

To illustrate the singularity structure of the relevant integrals, let us consider diagram (e). 
The corresponding amplitude squared has the form 

\M(q,k)\ 2 = 2C F C A g 4 - (2fc_ + g ) 

1 K Jl s k„q + + q_) (k + q_ + k„q + - 2kqcos8) v ' 

and the correction to the broadening is obtained by evaluating the integral 
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x 9(k + - k_) 6{q + - q_) e -*( k ±+i±>x± e -i(*+s) . (108) 

The exponentials arise from the Laplace-Fourier transform and the ^-functions ensure that 
both particles are emitted into the left hemisphere. Other than this, the two-particle phase 
space relevant for the soft sector is completely unrestricted, in contrast to the collinear phase 
space for the jet function, where the sum of the large momentum components had to add 
up to the center-of-mass energy Q. After integrating over the k + and q + components, we 
parameterize 

, riuvw riuvw ., riu riu /,™% 

k— — , q- = , \k± = — , \q± = — . (109) 

r r t t 

The integration over rj — . . . oo can immediately be done. Also the integration over w can 
be carried out. Due to the hemisphere constraint, this variable runs from i« = 0...-ifi><| 
and from w — ... e otherwise. Up to a prefactor, the integral then takes the form 

1 duu- 1 - a - 2e u- 1 - a - 2 * I 1 dv (v 2a 9(v - 1) + v 2a ea - v)) v l ~ a v- a - 1 {VU + 2U ^ 
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J v 2 + v 2 — 2vv cos 9 J \ J 

The terms in the first line have soft singularities for u = and it = 1 as well as for v — 1 (for 
other diagrams, also v = is singular). For v — |, the 6*- integral in the second line produces the 
collinear singularity shown in (I103p . As in the computation of the jet function, we introduce 
a parameter v cnt with < v cnt < \ to separate the soft and collinear regions. We then split 
the integration region into eight sectors as shown in Figure [9j In the four red corner areas, 
we disentangle the overlapping singularities in u and v using sector decomposition. In the 
blue collinear regions, we can either work with subtractions as in (11041) . or we can disentangle 
the singularities in the v and 9 integration using sector decomposition. We prefer the latter, 
since the subtraction technique becomes cumbersome for the self-energy diagrams (f) and (g), 
which involve the collinear denominator v 2 + v 2 — 2vv cos 9 to the second power. For the 
numerical evaluation, it is advantageous to combine the individual diagrams, since the sum of 
the diagrams is typically less singular than the individual ones. This is in particular the case 
for the self-energy diagrams. Working with the sum of the diagrams as given in Appendix [UJ 
we manage to perform the numerical evaluation with a relative accuracy of 10~ 6 or better. 

C Soft-gluon matrix element 

We write the squared matrix element describing the emission of two soft gluons in the form 
\M 2g (k,q)\ 2 = gt \c 2 A(k,q) + C F C A Af(k,q) 
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+ C F (T F n f - C A ) S 1 (k, q) + C F (T F n f - (1 - e)C A /2) S 2 (k, q) 



(111) 



The relevant diagrams are shown in Figure [3J The individual pieces of the amplitude squared 
are 
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The abelian part A(k, q) is just one half of the square of the one-loop amplitude. The non- 
abelian piece N{k,q) gets contributions from the box diagrams (a), (b) and (c), as well as 
from diagrams (d) and (e) involving a triple-gluon vertex. The latter involve a propagator 
denominator k ■ q. The last two structures S\(k,q) and S2{k,q) arise from the self-energy 
diagrams (f) and (g). 



D Results for individual diagrams 

In the main text, we have discussed the evaluation of diagram (b) in Figure [3] in detail. In 
the following, we summarize our results for the remaining diagrams. Using partial fractioning 
identities, we first relate diagram (a) to the other two box diagrams giving 
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and this relation carries over to the Laplace-Fourier transformed expression. 

We next consider diagram (c), which can be written as a convolution of two NLO integrals. 
In broadening-momentum space, the exact result reads 
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After Laplace and Fourier transformation this turns into 
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which is valid up to finite terms in the a-expansion. Here and below, we drop the index L on 
the variables z± = + z\ ± l)/4 for convenience. 



The calculation of the remaining diagrams proceeds along the same lines as in Section I4T2 
In broadening-momentum space, the non-abelian diagrams can be written in the form 
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and 
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7T 



,2c 



1+3/ £-l-2a 
d£ 

i-v 2-^ 



dp ^ 



-l-a 



(l + y-O-z-^-i + yY 



eV(l + p) + (2-0 2 (l + 2p) + 



i+p (i + P )[(2-o 2 + pe]-Vp 

e 2 (2-QV(i + P ) 

(2 - o 2 + P e 



(118) 



Note that the second integral generates a divergence in the analytic regulator in the limit 
p — > 0. Upon expansion in a, we could solve both of these integrals in terms of hypergeometric 
functions, but the exact expressions are rather lengthy In the limit y — > 1, we find 



i™(v i) ~ 2— r ( 2 - £ ) r (- £ ) + (i - „r' 2- 1 - r (i-2 £ )r(i + 2g ) 

h ( y ^ L >- Z r(3-2e) +( yJ r(l-e)r(l + e) 



x |l_l n (l-y)-ln2 + ^(^ + e) -*(l + e)J 



(119) 



and 



/ (2e) ( y 1) ~ 2— ara-On-O + (1 _ v w 2 -i-« r(i- 2e )r(i + 26) 
J5 r(i-2e) +(i y) 2 r(i-e)r(i + e) 



x |^-^ + ln(l-2/)-31n2-*^ + e) -3*(l + e)-4 Ts J 



(120) 
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This needs to be combined with the expansions for y < 1, which are given by 
lf d \y < 1) = ^ 2 ~ y arcsin y - 1 ln(l - y 2 ) - ^ In 



+ \{2U 2 ( l -^- y2 ) - 2U 2 (^ /l - y2 - 1 ) - 2 -LiJ- + ln 2 (l - y 2 

y 2 l 2 Vi + v /T— ^ 2 \i + y fr^2j y 2 l Vi + y^ 1 



+ ^ X 2y ^ U C h (tt - arcsin y) - Cl 2 (tt - 2 arcsin y) ^ - y 2 ln(l + ^1 - y 2 ) 

- 2 ln(l - y 2 ) ln(l + vT=7) + % ln(l - y 2 ) - In 2 ( 1 

2 8u V 



y vl + y 



3/ 2 + ln2 1 /l-y\ 2 tt 2 



2y VI + y/ 6y 



\n(—^\-y 2 -'-\e + 0{e 2 ) (121) 



and 



I? \V < 1) = -^j 1 + (W + \^?) - ln(l - */ 2 ))e + (^2Li 2 (i-^ 



i - r 



1,9/^ ■)> / _ On, /I ITn TT 



2 ln 2 (l + x/T^) + 2 ln2 (! - ^) - 2 ln(l - y 2 ) ln(l + V / W0 + y ) e 2 + 0(e 3 



- 2 In (-i+Vj-jg) - 1 ln(l - , 2 ) - A In (l±JL 
2 / y 2 % Vl -y 



+ 1{ 2 (1 - y 2 )Li 2 ( | V ] t \ - 2(1 + 2y 2 )Li 2 ( ) /l f 1 ) _ 2 y 2 l n2 ln (i _ ^ 



1 + 

- 6y Li 2 ( - - V ln 2 (l + v^^) - I In 2 In (1^) + ln 2 (l - y 2 



2(1 - y 2 ) ln(l - y 2 ) ln(l + - E±^±V ^ ^ 

8 V 1 + v 



+ % 2 In 2 ln(l + x/T^y 2 ) - |vr 2 je + 0(e 2 



(122) 



In Laplace-Fourier space the divergent terms in the ct-expansion become 

T (2d), . (OL s \ 2 C F C A( 2 -2\^{ U + f L) 2a 

x |- + -± + M2Li 2 (-— J + 81n 2 ^ + + 81nz + + 4 + -7r 2 J- + 16/i 2 (z L ) 

-8Si, 2 (- ^) + 4Li 3 (- ^ +8Si )2 (-tw) -40Li 3 (-w) -40S li2 (l - w) 
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+ 8Li 3 (l -w) +40S 1)2 
81n2 Li/ 



1 — w 



8 Lie 



1 — w 



( z +\ 



Z- \ 



^201n(l + 4) + 81n(4z + )jLi 2 ( 



-—) 

z - \ 

z + 



V64/ 
-w) + 8Li 2 



J + 3 lnZ+ 



+ 2 ln(l + z{) In 2 z + + 5 ln 2 (l + z l L ) In z+ + (16 - 28 In 2) In 2 z + + 8 In 2 ln(l + *£) In z + 



2 

+ 101n2 ln 2 (l + z\) + 8(2 + y -41n 2 2 - 2iu) In 2+ - 16(1 + 21n2)w 



+ 8(ln 2 2 + w) ln(l + 4) - 12 In 3 2 + 24 + tt 2 



3 3 s3 



and 



S L {T L ,T R ,Z L: Z R ) = y—j —^—(fj,T L ) {V+T L ) 

2 81n^+ 



2a 



(123) 



(124) 



"77 + 



8Li 2 ( ) +161n^+ + 7r 2 

z + 



a 1 



+ 



4 In z A 



(^28Li 2 {- Z —} + 16 In 2 z + + 27r 2 J i - 72 Si, 2 (- j-^j + 20Li 3 ( 



z_ 
1 — w 



+ 24Si >2 (-w) - 56Li 3 (-w) - 56S 1>2 (1 -w) + 24Li 3 (l - w) + 56Si >2 ( — ^) 



24LL 



1 — w 



( z + \ ( z -\ 



+ 72 In -r Li 2 



V 2 J 



z + 



j - 81n2 Li 2 



56 



In 3 z A 



24 In 2 Li 2 



1 — w 



28 ln(l + z 2 L ) + 24 In(4z+) U 2 (-w) + 6 ln(l + z\) In 2 2 



+ 71n 2 (l + z\)\xvz + - 52 In 2 In 2 z + + 24 In 2 ln(l + z£)ln* + + 14 In 2 ln 2 (l + ^ 

- f 64 In 2 2 - -tt 2 ) In z + + 24 In 2 2 ln(l + z\) - — In 3 2 + — In 2 - 6C3I - 1 , 
V 3/ 3 3 J^J 



where z± = (w ± l)/4, to = \Jl + z\ and the function h 2 (z) from (I6"4"j) . 

We next turn to the quark-loop contributions to the self-energy diagrams, which we write 
in the form 



(125) 



47T3-2E 



6? q U 



2a 



This time, we factored out the distribution (1 — y 2 ) 2 e . The two-dimensional integral repre- 
sentations of the diagrams read 



4 2/,n/) (y) = v e 



(1 _ y2) 3/2 ,1-H, 



7T 



2-f 



2-a 



1 + P 
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X 



2(1 - y 2 ) - £(2 - o (i + y - erHe - 1 + y)~* 



[(2 - o 2 + pe 2 ] [(i + pW - o 2 + pe\ - 4y* P ]< 



(126) 



and 



r( 2 3> n /) 



(y) = y 



2c 



(l_„2)3/a 



7T 



^ £l-2a ^ _ ^ 



X 



1 + 



e 2 (2-Q 2 (l+p) 2 

[(2 - o 2 + pe] 2 



2-{ 



2-Q 



(1+P) 2 

(i + y-o-5- e (e-i + y)-5- e 
[(i + P )[(2-o 2 +pe]-Vp] 2 - 



(127) 



As we factored out the distribution (1 — y 2 ) 2 e , we now have to expand to subleading order 
around the limit y — > 1. We find 



r (2/,"/) 



(!/->!) 



r(2-e)r(3-e) 
T(5 - 2e) 



and 



f 9 ' nf) (y 1) ~ 2- 2 " 2e r(2 ~ £) r(3 ~ £) 



T(5 - 2e) 



3e 
1 - e 



(i-y) 



(128) 



(129) 



For arbitrary y < 1 the expansions become 
1 - 3y 2 + 2y 4 



r(2/,n/) 



(y<i) 



■VT 



?p arcsm y 



(5-7y 2 )(l-y 2 )\ (1-y 



64y £ 



In 



1 + y 



+ 



3 - 6y 2 + 15y 4 - 16y 6 f 3 - % 2 + 37y 4 - 39y 6 , /l-y\ 3 - 6y 2 + 7y 



+ 



96y 4 

1 - 3y 2 + 2y 4 

8y 5 



+ 



192y £ 



In 



1 + y 



(5-V) ( i- y y/ 6 . 

128y 5 V 



3 - 6y 2 + 15y 4 - 16y 6 , 
H 3 hi 

48y 4 



1 — y 2 ( 4C12 (7r — arcsin y) — CI2 {tc — 2 arcsin y) 

/ 1 -y 



i-y 
i + y 



In" 



i + y 



+ 41n2 In 



48y 4 



i-y \ , 4 
i + y 



■a/i 3 ^ 



-7T 

3 



1 - y 2 \ 9 - 18y 2 - 51y 4 + 80/ 



288y 4 



e + 0(e 2 ) (130) 



and 



3(3 - y 2 ){l - y 2 f_ ^ (1 - y\ 27 - 54y 2 + 23y 4 



128y 5 



i + y 



192y 4 



27 - 54y 2 + 23y 4 , /1 + A/l-y 2 \ 99 - 233y 2 + 165y 4 - 39y 6 , /1-y 
hi ' 



96y 4 

,2\/i „,2\2 



1-y 2 



384y £ 



3(3- y -)(l-^)V 16Lii 
256y 5 V 



^ + In 2 f 



1 + 2/ 
39 



1 + y/ 



+ 4 In 2 In 



In . 

i + y 
i-y \ 

1 + y 



-7T 

3 



135 - 294y 2 + 13% 4 27 - 42y 2 + 19y A 



576y 4 96y 4 
In Laplace-Fourier space this finally translates into 



l-y^e + O^). (131) 



s'i""\r L , r R , z L , z R ) = (^)C F T F „, (ffff (132) 
4 /16 , 2 20\1 8, 2x , , . 8, , , 



16. / z_\ 32, o 80. 116 + 202? . 20 „ 472 2 



Lio In z + In z+ H w In z + z T ix 

3 \ z+J 3 + 9 + 9 + 9 L 27 3 

46 + 682? + IO2 4 , ,„ 2N 16 + 401n2 4 122 + 2721n2 2 110 + 1961n2] 

ln(l +zf)-\ Zj H zi H > 

9w v l; 9w l 9w l 9w J 

and 

S L {T L ,T R ,z Lt z R )= — C F T F n f (// t l ) (133) 

V47T/ aw; 

{-^-2(l + 4) 2 ln(l+ W )+ 1 + 6 '^ 

where the functions h\(z) and /i2(z) can be found in (1641 ). 

We proceed similarly for the gluon and ghost contributions to the self-energy diagrams, 
which we add up and write in the form 



S W/^\b L ,b R ,pi,Pi) (134) 

a 2 s C F C A ( /i 2 e 7fi ^ 26 1 fv+\ 2a n „,2\-2-e X d-2f„±\ X f h _\ T {2f/2g,C A ) l 

32vr 3 - 2e 

with 

(i-*/ 2 ) 3/2 f 1+y ^ ,1- 

'l-y 



lf f ' CA \y) = 4y 2 * (1 ~ ^ f +V di e~ 2a (2 - (1 + 2/ - 0"^(£ - 1 + y)^ (135) 

^ Jl-y 



f^H P l ~ a ( 2 ~ ^ + ~ Q P( 2 ~ + I 2 / 3 - + £ (1 ~ P)[( 2 ~ 2 ~ ^ 

io 9 1 + P [(2 - 2 + P£ 2 ] [(1 + p) [(2 - 2 + PC 2 ] - 4y 2 p] 2 

and 

i ( s 9 ' CA \y) = y 2 ' ^~ rr ' I * ' "de ^ 2Q (2 - f^d, 



tt A-/^ ( ' 4 U p [(i + P)[(2-e) 2 + pe 2 ]-4y 2 P p 



x p^ Q 



(l + e )(l + p 2 ) + (6-2 e )p , (l + e )[(2-0 4 + p 2 e 4 ] + (6-26) P e 2 (2-0 



(1 + P) 2 [(2-0 2 + P^ 



212 



(136) 
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The solutions in the limit y — > 1 now read 



/f^(^l)^-2-- r(1 - e)r(3 - e 



-2- 1 -* (1-j/) 



r(5 - 2e) 

T(l - 2e) r(l + 2c) 



5-3e- (43-21e-4e 2 ) (1 - y) 



r(2-e)r(l + e) 



(137) 



and 



T(5 - 2e) 



f^)(^l).2— ^| )r{3 - 6 

-2-e M M-e r(l - 2e) T(l + 2e) 

r(2-e)r(i + 6) 



5 - 3e - (3 - 5e + 4e 2 ) (1 - y) 



+ 2-^ (1 - y) ] 



(138) 



and the expansions for y < 1 are 

< 1) ^l-VKl-^^.^ + (5 -16^(1-^1 - , 



2» 5 



3 - 9t/ 2 + 36t/ 4 - my'' < 3 - 9y 2 + 16y 



2V 

(1-V)(1-^ / 



16y 5 VI + y 

3 - 3y 2 - 22y 4 + 27y 6 



12|/ 4 

— arcsin y) — CI2 (fl" — 2 arcsin y) 



In 



1-2/ 



2y 5 



12y 5 

(1-^)3/2 



2/ 



In 2 



1 + y 



41n2 In 



2y 5 

l + y 



arcsin y 
4 

-7T' 

3 



+ 



3 - V + 36y 4 - 40y b , /1+Wl-y 



I2y* 



In 



9 - 18y 2 - 84y 4 + 124y € 
3%4 



e + C(e 2 ) (139) 



and 



ip' CA \ y < i) = (9-V)a-V) 2 i n ( l - y\ ( 9 - 4 ^ 3 - 5 ^ 2 ) 



32y 5 



1 + 2/ 



48y 4 



(9-4y 2 )(3-5y 2 ) , (\ + v 7 ! ~2/ 2 \ 18 - 44y 2 + 33y 4 - 12y 6 , /l 
In ' 



24y 4 

(9-4y 2 )(l-y 2 ) 2 / 
64?/ 5 



16Li 2 



24j/ 5 



In 



1 + 2/ 



1-2/ 
1 + 2/ 



In 2 



1-2/ 
1 + 2/ 



41n2 In 



1-2A , 4 
1 + 2/ 



-7T 

3 



27- 66y 2 + 8y 4 27 - 45y 2 + 28y 4 f -\ 2 

727 24? VW) C + 0(0- 



(140) 



In Laplace-Fourier space we finally obtain 



t(2/,CU). 



S L {t l ,t Ri z l ,z r ) = i^— 



i - \ 2a 



a 



(141) 
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f 10 /40 , 5 62\1 8, 2N , . . 8, 2x , , , 5 , 

40 T / z\ 80, 2 248, 278 + 20z? , 266 + 5961n2 2 
+ — Li 2 1 I + — In z + + — In z + w In z + z L 

3 V z+/ 3 9 9 9w 

109 + 149z? + lOsf , ,„ 2 , 16 + 401n2 4 266 + 4661n2 20 2 12221 

H ln(l + Zj) z A T 1 z r H \ 

9w 9w L 9w 9 27 J 



and 

5 L (ti, 7*, ^) = (- j — (a* rj - 

x || + (1 + 4) (3 + 2,2) ln(l + W ) + 1 " lhZ \ ~ 64 ln(l + z\) 

+ ^-51n2 + ^- + 2^) W |. (142) 
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